function res = gauss_lin(p,x)
   res = p(1)*exp(-(x - p(2)).^2 ./ (2*p(3)^2)) +  x*p(4) +p(5) ;
    
end